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Magnetic field amplification in CDM anisotropic collapse 
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ABSTRACT 

We use the Zel'dovich approximation to analyse the amplification of magnetic fields 
in gravitational collapse of cold dark matter during the mildly nonlinear regime, and 
identify two key features. First, the anisotropy of collapse effectively eliminates one 
of the magnetic components, confining the field in the pancake plane. Second, in 
agreement with recent numerical simulations, we find that the shear anisotropy can 
amplify the magnetic field well beyond the isotropic case. Our results suggest that the 
magnetic strengths observed in spiral and disk galaxies today might have originated 
from seeds considerably weaker than previous estimates. 
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1 INTRODUCTION 

Large scale magnetic fields, with strengths 10 _7 -10 -5 G, 
have been observed in spiral and disc galaxies, galaxy clus- 
ters and high redshift condensations [see (Kronberg 1994; 
Han & Wielebinski 2002) and references therein]. The most 
promising explanation for the large scale galactic fields has 
been the dynamo mechanism, with the required seeds com- 
ing either from local astrophysical processes, such as bat- 
tery effects, or from primordial magnetogenesis [for recent 
reviews, see (Grasso & Rubinstein 2001; Widrow 2002)]. The 
linear evolution of large scale magnetic fields and their impli- 
cations for structure formation has been studied by several 
authors [see, e.g., (Ruzmaikina & Ruzmaikin 1971; Wasser- 
man 1978; Papadopoulos & Esposito 1982; Evans & Fennelly 
1985; Tsagas & Barrow 1997; Barrow & Subramanian 1998; 
Tsagas & Maartens 2000)]. Certain aspects of the mildly 
nonlinear clustering can be analysed in spherical symmetry, 
but this approximation inevitably breaks down as the col- 
lapse proceeds and any initially small anisotropics take over. 
When magnetic fields are involved, the need to incorporate 
anisotropic effects is particularly important, as the fields are 
themselves generically anisotropic sources. Numerical simu- 
lations of magnetic field evolution in galaxy clusters suggest 
that anisotropies lead to additional amplification of the field. 
Tidal effects during mergers, for example, increase the mag- 
nitude of the magnetic field (Roettiger et al 1999). Shear 
flows in galaxy clusters can also amplify the field beyond 
the limits of spherical compression (Dolag et al 1999; Dolag 
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et al 2002). The latter simulations also suggest that the fi- 
nal magnetic configuration is effectively independent of the 
field's initial set up or of the presence of a cosmological con- 
stant. 

Here we use the Zel'dovich approximation to look an- 
alytically into the effect of gravitational anisotropy on 
seed magnetic fields in the mildly nonlinear regime. The 
anisotropic collapse is driven by cold dark matter (CDM). A 
previous analysis (Zel'dovich et al 1983), where the matter 
is purely baryonic, concluded that the seeds must be negli- 
gibly small in order to avoid the growth of magnetic fields 
to levels which prevent pancake formation. By contrast, this 
strict constraint is avoided when CDM dominates, since the 
magnetic field couples only gravitationally with the CDM. 
The field is frozen into the baryon fluid, and baryons are 
dragged by the CDM gravitational field. The baryons are 
in a different state of motion to the CDM and, unlike the 
CDM, they feel the magnetic backreaction. As a first ap- 
proximation, however, we will ignore the magnetic backre- 
action and the relative motion of the baryons, effectively 
considering a single fluid. This approximation also maintains 
the acceleration-free and irrotational nature of the motion, 
which are key ingredients of the Zel'dovich approximation. 
Our approach may be seen as a qualitative starting point 
for a more detailed analysis. 

We begin in Section 2 with a dynamical system descrip- 
tion of the Zel'dovich approximation, which directly shows 
that, in a generic collapse, pancakes are the (local) attrac- 
tors (Bruni 1996). In Section 3 we consider the dynamics 
of the magnetic field as it collapses with the matter. The 
magnetized dynamical system is five-dimensional. Pancakes 
are still the attractors, with the magnetic field squeezed in 
the pancake plane. Note that, as the galaxy is formed, tidal 
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forces are generally expected to change the orientation of 
the galactic plane relative to the pancake. Nevertheless, the 
confinement of the field in the pancake plane that we find 
here, is qualitatively consistent with magnetic field observa- 
tions in numerous spiral and disc galaxies. We also provide 
some quantitative results by relating the growth of the field 
to that of the matter density contrast. When comparing the 
shearing to the isotropic collapse, we find that anisotropy 
can lead to an appreciable increase in the amplification of 
the magnetic field. These analytical results are in qualita- 
tive agreement with those of the earlier mentioned numer- 
ical simulations. We interpret them as an indication that 
the magnetic strengths observed in numerous galaxies and 
galaxy clusters today could have resulted from seeds consid- 
erably weaker than previous estimates. 



cannot consider a purely initial value problem. Instead, as 
a consequence of action at a distance, one necessarily needs 
boundary conditions (Barrow & Gotz 1989). 
The Zel'dovich ansatz, 



-V x <p 



e 



V72 - 



-Ei 



(7) 



implies that the parentheses in (5) and (6) vanish, leading to 
a local system of ordinary differential equations by eliminat- 
ing the dependence on S and Eij. Since the shear and tide 
matrices now commute, the above system can be written in 
the shear-tide eigenframe and reduces to three equations, 
one for 6 and two for the independent shear eigenvalues ai , 
(T 2 . It is then straightforward to verify that the solutions of 
the reduced equations that follow from the Zel'dovich ansatz 
are 



2 ZEL'DOVICH DYNAMICS AS A 2-D 
SYSTEM 

The Zel'dovich approximation (Zel'dovich 1970) [see 
also (Shandarin & Zel'dovich 1989; Padmanabhan 1993) for 
a review and (Buchert 1996; Matarrese 1996a; Matarrese 
1996b; Ellis & Tsagas 2002) for an approach similar to ours] 
arises from a simple ansatz, which extrapolates to the non- 
linear regime a well known result of the linear perturbation 
theory. In the Eulerian frame (x,t), which is comoving with 
the background expansion, one defines the rescaled peculiar 
velocity field u = dx/da, where a — a(t) is the background 
scale factor. Then, to linear order and ignoring decaying 
modes, u = — V x <^, where (p = 2y>/3-ffo a o i s the rescaled pe- 
culiar gravitational potential and H = a/a. The Zel'dovich 
approximation uses this linear result in the rescaled nonlin- 
ear continuity, Euler and Poisson equations, which are 



85_ 
da 
du 

da 



= -Vx • (<5u) - VxG , 

3 3 

= -u- V x u- — u — V x <£, 

2a 2a 

= u, 



(1) 

(2) 
(3) 



where 5 — Sp/ph is the density contrast. Assuming irro- 
tational motion, we define the peculiar expansion and the 
peculiar shear by 6 = V x • u and Oij = 
spectively (di = d/dxi) 
Eqs. (1) and (2) give 



Then, using convective derivatives, 
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- d'ikS'kj ~ —(d'ij + Eij) , (6) 



djdi(p — (d 2 p/3)Sij is the Newtonian traceless 
tidal field and 2a — aijCTij. One can then substitute d ip 
from (3) into (5). In the linear, or in the shear-free case, 
Eqs. (4) and (5) form a local system describing the evolu- 
tion of the fluid along its flow lines. In the general case, the 
presence of the tidal term Eij in (6) and the lack of a cor- 
responding evolution equation show that the above system 
is non-local and not closed. It also emphasises the fact that 
in Newtonian gravity, as opposed to general relativity, one 
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where Ai are the three eigenvalues of the initial tidal 
field (Matarrese 1996a). In particular, <5 zcl is the solution 
of the continuity equation (4), when O is given by (8). 
On the other hand, if <5 dyn is the density contrast used 
in the Poisson equation, then from (3) and (7) one gets 
<5 dyn = — a6 zo1 / <5 zcl as a consequence of the approxima- 
tion (Nusser et al 1991). Note that a negative eigenvalue 
Ai corresponds to collapse along the associated shear eigen- 
direction. Also, 1-dimensional planar pancakes are solutions 
corresponding to two vanishing eigenvalues. For example, 
Ai = = A2 and A3 < describes 1-dimensional collapse in 
the third eigen-direction. As is well known, this is not only 
a solution of the simplified local dynamics that follows from 
the Zel'dovich ansatz, but also an exact solution of the full 
system (3)-(6). In general one expects at least one of the Ai 
to be negative and smaller than the other two, and on this 
basis the generic solution should to tend to a pancake. 

The existence of pancake attractors can be confirmed 
by a dynamical system approach (Bruni 1996). Dropping 
the superscript "zel" , we define the new time variable 



/ 



Oda , 



and the dimensionless dynamical variables 



E+ = _(Ei + E 2 ) : 



V3 



(Ei-E 2 ), 



where Ei = di/®, to arrive at the system 



(11) 
(12) 
(13) 

3 L- -v-..-, .- TJ — (14) 

Equation (5) has now decoupled and the problem has been 
reduced to the study of the planar dynamical system (13), 
(14) depicted in Fig. 1. The evolution of the shape of a 
generic fluid element is given by the time dependence of its 
axes £j. In the shear eigenframe (Ellis 1971), 
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Table 1. The fixed points in the £+, £_ plane. The second and 
third columns give a typical example in the case Si = £ 2 on 
the line £_ = 0, the fourth and fifth give the eigenvalues of the 
Jacobian of the dynamical system (13), (14) at that fixed point. 



Fixed Point 


£+ 


£_ 




s_ 


Stability 


pancakes 


-1 





-1 


-1 


asymptotically 
stable node 


filaments 


1 

2 





1 

2 


1 

2 


saddle 


spherical 








1 

3 


1 

3 


unstable node 



dli_ 
da 



(I 



(15) 



One-dimensional pancakes are locally axisymmetric solu- 
tions of the planar system with two of the Ei equal to —1/3, 
so that two of the are constants. Generic pancakes on the 
other hand, correspond to solutions that are asymptotic to 
one of the three 1-dimensional pancakes (one for each shear 
eigen-direction) . In this case two of the £; tend to constants 
(i.e. two of the Ei approach —1/3). Exact 2-dimensional "fil- 
aments" are solutions with two Ei equal to 1/6 (e.g. E_ = 
and E+ = 1/2). 

The system (13), (14) admits two sets of invariant sub- 
manifolds, corresponding to straight lines in the E + ,E_ 
plane of Fig. 1. The two sets of lines in Fig. 1 are: (i) those 
that form the central triangle (e.g. E3 = —1/3, E+ = 1/2), 
representing trajectories with one of the Ei = —1/3; (ii) 
those that bisect the vertex angles (e.g. the E + axis, E_ = 0) 
represent locally axisymmetric solutions with two equal Ei. 
The vertices of the triangle are the stationary points of 
(13), (14) that correspond to exact 1-dimensional pancakes. 
Where the first set of lines intersects the second, we have 
stationary points representing exact filamentary solutions 
(e.g. E_ = 0, E+ = 1/2). Finally the centre, where the three 
bisecting lines meet, corresponds to spherical collapse. It is 
clear from Fig. 1 that depending on the initial conditions, a 
generic fluid element collapses to one of the three pancakes, 
while the three filaments are unstable. Indeed, looking at 
the eigenvalues s± of the Jacobian of (13), (14), one finds 
that the pancakes are asymptotically stable nodes, while the 
filaments are saddle points, and the spherical solution is an 
unstable node (see Table 1). This proves that, once collapse 
sets in, pancakes are the attractors of the generic trajecto- 
ries. 



3 SHEAR AMPLIFIED MAGNETIC FIELDS 

We now analyse the effects of mildly nonlinear clustering 
on magnetic fields frozen into a collapsing protocloud that 
is falling into a CDM potential well. We neglect magnetic 
backreaction on the fluid (effectively we consider a force- 
free field), as well as the velocity of the baryons relative to 
the CDM. In comoving coordinates the magnetic induction 
equation reads 
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Figure 1. Phase plane £4. = X, £_ = Y. The lines forming the 
central triangle correspond to one of the Ej = —1/3. The three 
pancakes are located at the intersections of these lines. The lines 
bisecting the vertices correspond to two equal £j. 
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Introducing the rescaled magnetic field vector Bt — a 2 Bi, 
and using the shear eigenframe, this gives 



dgi 
dr 

di?2 
dr 



-i(E+ + ^E_-2) B 1: 
-i(E+-v/3E_-2) B 2 . 



(17) 
(18) 
(19) 



A pancake solution, with collapse along the third shear 
eigen-direction, is characterised by the asymptotic values 
E + = — 1 and E_ = 0, so that 

£>i,2 oc a~ 2 e T oc a~ 2 exp(— J Qda), B3 oc a~ 2 . (20) 

In the above a~ 2 gives the dilution due to the expansion, e T 
the increase caused by the collapse and r provides a mea- 
sure of the collapse timescale. Thus, we expect the magnetic 
strength to increase as the field collapses with the fluid and 
gets confined in the pancake plane. This qualitative result is 
consistent with the pattern of the magnetic fields observed 
in spiral and disc galaxies. 

To obtain a more quantitative result we need to solve 
the magnetic evolution equations (17)-(19). Using expres- 
sions (8) and (9) we arrive at 
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'(l+aoAi)(l + a A 2 ) 


Bs 


(l + oAi)(l + oA 2 ) 



00 y 

a / 



/ao 
V a 



(21) 
(22) 
(23) 



where ao, Bi a are initial values. In 1-dimensional planar 
collapse along the 3rd eigen-direction (Ai = = A2 and 
A3 < 0), the pancake singularity is reached as a — > — I/A3. 
In that case B3 decays, whereas B\, B 2 increase without 
bound (neglecting backreaction). 
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We now compare with spherical collapse, limiting to the 
Zel'dovich case which is given by three equal eigenvalues. 
Then, any one of (21)-(23) reduces to 



(24) 



B sph = / l + a t ,A \ 2 /oo\ 2 
Bo \ 1 + a\ J \a J ' 
and all the magnetic components diverge as we approach 
the point singularity, a — > — 1/A. For comparison, we assume 
equal initial values Bi = Bo, which leads to the ratio 



gi 

iJsph 



(l + a„A2)(l + aoA 3 )(l + aA) 2 



(l + aA 2 )(l+aA 3 )(l+a A) 2 ' (25) 

with an analogous relation for B2/B Bph . This result shows 
that the relative growth of the field depends on the eigenval- 
ues corresponding to the two alternative patterns of collapse. 
When A3 < A < 0, or equivalently when < — I/A3 < —1/A, 
the anisotropy dominates the collapse. Thus a — > — I/A3 be- 
fore a -> -1/A, and therefore Bi, B 2 > B sph . In other 
words, anisotropy can lead to a considerable increase in the 
strength of the magnetic field. 

Conditions such as A3 < A arise naturally whenever 
we compare pancake to spherical collapse. Indeed, consider 
pure pancake collapse along the 3rd shear eigen-direction 
(Ai = = A2). Then, from (8) we have 



e 



a 3 = 

1 — a 9o 

For spherical collapse on the other hand, Eq. (8) gives 

e 



(26) 



(27) 



3 — a 6o 

so that A3 < A. Accordingly, one should expect stronger am- 
plification for magnetic fields frozen into anisotropic, rather 
than isotropic, collapse. Note that we have implicitly as- 
sumed that, initially, the value of is the same both for the 
pancake and for the spherical collapse. 

An alternative way of estimating the relative amplifi- 
cation is via the initial density contrast. For Ai = = A2, 
Eq. (25) leads to 

(l + a A3)(l + aA) 2 



Bi 

B s P h (l + aA 3 )(l + a A) 2 



(28) 



with the same relation for B2/B sph . Substituting for A, A3 
from (26) and (27) we arrive at the expression 

Bi,2 [1 - <W(1 + &,)][! -aWa (3 + fr)] 2 . , 

B** h [l~a8 /ao(l + So)][l-So/(3 + 5 W ~ ' V ' 
where the equality on the right hand side holds at a = a . 
The above directly relates the ratio Bi^/B sph to the over- 
density So at the onset of the nonlinear regime. In deriving 
Eq. (29) we have also used the relation 80 = — aoOo, which 
arises directly from the Zeldovich ansatz, Eq. (7). Accord- 
ing to (29), the ratio Bi/B sph keeps increasing as the col- 
lapse proceeds towards the pancake singularity, namely as 
a -> o (l + <5,7 1 ). At that point Bi, 2 /B sph diverges. Note 
that, for the same initial O, the isotropic collapse reaches 
the singularity at a — a (3 + S^ 1 ), that is later than the 
pancake collapse. On using Eq. (29) we can calculate the 
relative amplification of the magnetic field, for a given 80, 
at any time prior to the pancake singularity. To obtain a 
numerical estimate consider the following example. Assume 
that 80 = 1 at the onset of the nonlinear regime. This means 




Figure 2. Typical growth of the ratio R = B li2 /B sph , Eq. (29), 
during the late stages of the collapse. The horizontal axis mea- 
sures the percentage of the total pancake collapse according to 
(a - a )/a X 100% (a < a < 2a for 8 = 1). 

that the pancake singularity is reached at a — 2ao, whereas 
the point singularity at a = 4ao. Then, the anisotropically 
collapsing magnetic field is one order of magnitude stronger 
than an isotropically collapsing one (i.e. Bi j2 — 10B sph ) by 
the time a — 1.95ao (i.e. when we have undergone 95% of 
the total pancake collapse). This means that the effect of 
the anisotropy on the amplification, though weak early on, 
increases fast as we approach the final stages of the collapse 
(see Fig. 2). Then, the presence of shear can strengthen the 
magnetic field well beyond the limits of simple isotropic com- 
pression. In other words, as the field is dragged deep into the 
potential wells of the CDM, one should expect an apprecia- 
ble increase in the magnetic strength due to shearing effects 
alone. 

Although Eq. (29) has been derived in the idealized case 
of a 1-dimensional pancake, with the spherical case treated 
via the Zel'dovich approximation, we expect it to be approx- 
imately true in more realistic situations as well. This rea- 
sonable expectation is further strengthened by the fact that 
our analytical estimates are in good agreement with those 
obtained by numerical simulations [see (Dolag et al 1999)]. 
These simulations also argue for an extra order of magnitude 
increase in the strength of the magnetic seed field, during 
the final stages of the collapse, because of shearing effects. 

Finally, we can gain further insight into the relative 
amplification of anistropically versus isotropically collapsing 
magnetic fields by rewriting (28) as 



gi 



ii f t ph \ 2 _ ti /2 £3 y 



(30) 



where we have also used Eqs. (15), (26), (27), and linear ini- 
tial conditions for two initially spherical fluid elements. This 
relation shows that the ratio Bi /B sp grows almost linearly 
with the eccentricity of the pancake I1/I3, independent of 
initial conditions. 



4 CONCLUSIONS 

Protogalactic clouds do not collapse isotropically: in any re- 
alistic situation one expects small anisotropies in the initial 
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velocity distribution, which will be amplified as CDM col- 
lapse progresses. We used the Zel'dovich approximation to 
investigate the effects of such gravitational anisotropy on 
the evolution of a seed magnetic field frozen into the baryon 
fluid. We considered the mildly nonlinear regime and ig- 
nored the magnetic backreaction on the baryons as well as 
the relative velocity between baryons and CDM. The CDM 
dominates the collapse and the gravitational anisotropy that 
amplifies the magnetic field. Our qualitative analysis shows 
that, as a generic result of the anisotropy of the collapse, 
the field effectively "loses" one of its components and is con- 
fined in the plane of the pancake. Although tidal forces are 
generally expected to change the orientation of the galac- 
tic plane relative to the pancake, our qualitative picture is 
consistent with magnetic field observations in numerous spi- 
ral and disk galaxies. More quantitatively, our results show 
a much more efficient magnetic amplification in a shearing 
rather than in a shear-free collapse. This analytical result 
agrees with numerical simulations of shear and tidal effects 
on the evolution of magnetic fields in galaxies and galaxy 
clusters. Taken together, these results indicate that, when 
the anisotropic effects of CDM collapse are taken into ac- 
count, the magnetic strengths observed in galaxies today 
could have originated from seeds considerably weaker than 
previous estimates. 
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